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The prior knowledge of the gravitational waveform from compact binary systems makes matched 
filtering an attractive detection strategy. This detection method involves the filtering of the detector 
output with a set of theoretical waveforms or templates. One of the most important factors in 
this strategy is knowing how many templates are needed in order to reduce the loss of possible 
signals. In this study we calculate the number of templates and computational power needed for 
a one-step search for gravitational waves from inspiralling binary systems. We build on previous 
fSJ ■ works by firstly expanding the post-Newtonian waveforms to 2.5-PN order and secondly, for the 

' first time, calculating the number of templates needed when using P-approximant waveforms. The 

, analysis is carried out for the four main first-generation interferometers, LIGO, GEO600, VIRGO 

(N . and TAMA. As well as template number, we also calculate the computational cost of generating 

5h ' banks of templates for filtering GW data. We carry out the calculations for two initial conditions. 

In the first case we assume a minimum individual mass of 1 Mq and in the second, we assume a 
minimum individual mass of 5 Mq. We find that, in general, we need more P-approximant templates 
to carry out a search than if we use standard PN templates. This increase varies according to the 
order of PN-approximation, but can be as high as a factor of 3 and is explained by the smaller 
span of the P-approximant templates as we go to higher masses. The promising outcome is that 
for 2-PN templates the increase is small and is outweighed by the known robustness of the 2-PN 
■ P-approximant templates. 
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The inspiral phases of close binary systems are expected to be the most important sources of gravitational waves 
(GW) for the first generation of detectors, LIGO, GEO600, VIRGO and TAMA Under radiation reaction the 

orbit of the binary decays emitting a signal which increases in amplitude and frequency called a "chirp" signal. The 
three main types of compact sources of GW interesting for detection purposes are neutron star-neutron star (NS-NS), 
' neutron star-black hole (NS-BH) and black hole-black hole (BH-BH) binaries. To see NS-NS binaries it is estimated 
rr' that we need to look to 100-400 Mpc to see three events per year While it is assumed that there is a greater 

population of NS-NS binaries, it is the BH-BH binaries which are the strongest candidates for detection from a purely 
bandwidth point of view. It is believed that with an initial network of detectors, it will be possible to detect 2^3 
events per year at a distance of 200 Mpc Q. Due to the masses involved, the inspiral phase for these systems shuts 
■ off earlier than in BH-NS or NS-NS systems. This is advantageous as these binaries should "chirp" in or about the 
most sensitive frequency of the detectors. 

As we have some prior knowledge on what the waveform should look like, the optimal strategy is "Wiener" or 
"matched" filtering j|] . This works as follows : we firstly create a set of theoretical models or "templates" based on a 
certain parameter set. These templates are then cross-correlated with the detector output. The signal-to-noise ratio 
(SNR) is compared to a predetermined threshold to decide whether or not a signal is present. The cross-correlation 
is weighted with the inverse of the noise power-spectral-density of the detector. This has the effect of emphasising 
the most sensitive frequencies of the detector. This is another reason why the BH-BH binaries will indeed be the best 
candidates. 

In order that the search for GW is not computationally exhaustive we need to have some predetermined criteria for 
the generation and placement of templates. In reality a template is dependent on a number of parameters such as mass, 
spin, eccentricity etc. of the binary. As the composition of the signal is unknown, the templates have to be laid out in 
a lattice such that the correlation between a possible signal and a nearby template achieves a predetermined threshold. 
In this study, the templates used represent non-spinning circular binaries. We therefore have a reduced parameter 
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space as we can theoretically model our system by using only the individual masses of the binary components. It has 
been shown that for detection purposes it is sufficient to work with restricted post-Newtonian waveforms |llo| . This is 
when the amplitude of the waveform is kept at the dominant Newtonian level, but the phase is extended to various 
post-Newtonian (PN) orders |ll|-|2l[]. While it may still be sufficient to work with restricted waveforms, there may be 
evidence to suggest that the standard PN waveforms may not be good enough. 

The PN approximation is a perturbative method which expands the equations of motion, binding energy and GW 
flux as a power series in v/c, where v is the orbital velocity and c is the speed of light. In the early stages of an 
inspiral, the radiation reaction time-scale, Trr is much greater than the orbital time-scale, Torb- It is in this adiabatic 
regime that the PN approximation works best. The main problem for PN waveforms comes when we enter the strong 
field regime. To get an idea of how well the PN waveforms perform, we can look at the ratio Trr/Torb for various 
systems in and around the most sensitive frequency of the detector. If we take the LIGO noise curve, we are mostly 
interested in the region between 100 — 300 Hz. For a 1.4-1.4 Mq binary the ratio between the time-scales lies between 
4.5 X 10'^ < Trr/Torb < 680. Howcvcr, for a 10-10 Mq binary, the ratio lies between 140 < Trr/Torb < 40, where the final 
value is taken to be the frequency at the last stable orbit {fuo ~ 190 Hz). We know that the adiabatic approximation 
breaks down when the two time-scales are comparable. We can see that even for the lowest mass BH-BH binaries 
expected to be detected, we cannot have much confidence in theperformance of the PN templates. Recently, Damour, 
Iyer and Sathyaprakash introduced P-approximant waveforms P3| . These templates are constructed by taking the PN 
expansion, defining new energy and flux functions and resumming the series expansion of the new functions using the 
method of Pade approximation. These new waveforms were shown to have better behaved convergence properties than 
the PN templates. This eliminated the problem with PN waveforms that a higher approximation didn't necessarily 
mean a better template. 

For the placement of templates it has been shown that by incorporating techniques of differential geometry we 
can greatly simplify the problem p^ . As the matched filtering is to be carried out in the Fourier domain, it can be 
shown that the restricted PN Fourier transform waveforms (RPNFT) can be thought of as vectors in a Hilbert space. 
Sathyaprakash has shown that when considering two-mass parameter waveforrns, a convenient coordinate system 
to use is not the individual masses but the Newtonian and 1-PN "chirptimes" |2^. The chirptimes are convenient 
combinations of the total mass, m = mi + m2, and the reduced mass ratio, 77 = mim2/m^ , of the system. The 
concept of a metric tensor on the search manifold was defined by means of two different but ultimately equivalent 
methods |2^,|2^ . The number of templates needed for a one-step search has already been calculated at Newtonian , 
1-PN and 2-PN ||2^ orders. In this paper we will extend the previous works in three ways : (i) use the more 
recent models of the noise power spectral density of the detector |^ (ii) expand the template number calculation to 
2.5-PN for RPNFT waveforms and (iii) calculate the template number for P-approximant waveforms to 2.5-PN order. 

The rest of this paper is organized as follows. In Sec. O we recap and partially extend the geometric method of 



GW detection. In Sec. [II we define the T and P-approximant waveforms, while Sec. IV contains a description of 
the numerical method used. Sec. ^ contains the results for the template number and computational cost using both 
types of waveforms for all ground-based detectors. 

Throughout this paper we use geometric units (G — c = 1). 



II. THE GEOMETRIC METHOD 



The restricted PN time domain waveform can be written as 

/i (t. A'') = a (i) cos $ (t, A^) . (1) 

where {A''} is an open set on M". As the search for GW will be carried out in the Fourier domain, it is necessary 
to have such a representation of the waveform. In geometric terms the Fourier transform is a bijection from the real 
space M" to the complex space C". As the map h{t, X'^) 1-^ h{f, A^) is diffeomorphic, this defines a topology on both 
manifolds [^0| . The phase of the GW in the Fourier domain can be conveniently parameterized by the parameter 
set A'^ = {tQ,^o,Tk}, where to is the time-of- arrival of the wave, $0 is the phase of the wave at time-of-arrival. The 
"chirptimes", {rfc}, are functions of the individual masses of the system defined to 1.5-PN approximation by 
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FIG. 1. The one-sided effective noise spectrum for the four ground-based interferometers. 



The parameter /o corresponds to the lower frequency cutoff of the detector and each subscript, k, corresponds to a 
(fc/2)-PN approximation. As we could, in principle, have used the parameter set {to, $0j ™ii 1^2} to represent the GW 
phase, we can see than only two of the chirptimes can be independent. This will prove convenient as it will allow us 
to choose a combination of any two chirptimes to serve as the coordinate basis for our waveforms. As the chirptimes 
are dynamical parameters dependent on the astrophysical system, we will refer to them as being "intrinsic" . We will 
thus refer to the kinematical parameters tg and $0 as "extrinsic" parameters. 

By treating the Fourier domain waveforms as vectors in a Hilbert space, we can define a scalar product 

{a\b)^2 ^^~a{f)b*{f) + c.c., (5) 

where c.c. denotes complex conjugate, Sh{f) is the one-sided noise power spectral density (PSD) and d{f) is the 
Fourier transform of a{t) defined by 

/oo 
dta{t)e'^'fK (6) 
-00 

We can provide an analytical estimate for the shape of the PSD. In Fig|l| we have plotted the effective noise \/ / Sh{f) 



for the four initial ground-based detectors |29|. For our purposes, we can think of the factor 1/ Sh{f) as simply a 



weight in the scalar product. We can also define a vector norm 

\x\^{x\x)"\ (7) 

It is possible to define a metric tensor on the manifold by using two separate, but ultimately equivalent methods. For 
both methods we assume that each waveform is normalized to unity. We also consider two points on the manifold 
separated by the infinitesimal coordinate distance AA^ In the first method we use the fact that the proper 

distance between the two points is given by 

ds^^\h{f,X>^ + ^\>^)-h{f,\>^)\\ (8) 



This can be rewritten as 



where the second expression comes from Eqn (|^) after we neglect all terms of order AA^ and higher. The quadratic 
form on the right hand side suggests we define the metric as 
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dh \ „ /■°° df dh* dh 



To fully justify this, we can use the fact that in any normed Euclidean vector space, the right hand side of Eqn (g) 
defines a metric pO|-^ . 

In the second method we start by calculating the scalar product between two templates h{\^) and h{X^ + AX^) 
separated by a distance AA^ in the parameter space. The scalar product is nothing more than the cosine of the angle 
between the waveforms. It has been shown that by expanding the scalar product to quadratic order we obtain p4[ 



dh 



AA'^AA^ (11) 



= 1 - -e^, AA^AA^ (12) 



The quantity on the left hand side above is referred to either as the match or overlap between templates. For this 
study we will use the term overlap. We thus rewrite the preceding equation as 

= l-^ds\ (13) 

where O denotes the overlap. As all the parameters of the waveform are not independent, the search for GW will 
only need to be carried out on the 2-D sub-manifold of independent intrinsic parameters. In order to project from our 
initial n-D space to the 2-D sub-space we recursively use a projection onto orthogonal sub-spaces. This is achieved 
by using 

7m.=?^.-%^. (14) 

This projection operation is equivalent to a mapping from the n-D space onto an orthogonal hypersurface. The 
projection itself has two separate meanings. From a geometric point of view, it can be seen as a minimization of 
distance between two templates. From a data analysis point of view, it can be viewed upon as a maximization over 
extrinsic parameters. By using this projection, we ensure that the overlap between templates on the 2-D sub-space 
is maximized over all extrinsic parameters. In Appendix A we outline some of the properties of this projection from 
both a geometric and detection point of view. 

Once we are on the sub-manifold of intrinsic parameters, the number of templates is given by simply dividing the 
proper volume of the manifold by the volume of the template, dP, i.e. 

where g — det\gij \ is the determinant of the metric tensor on the 2-D manifold. For a rectangular lattice of templates 
in N dimensions, the volume of each template is given by 

dl^= ^2^^^-1111 j . (16) 

where MA4 is the minimal match. This is a quantity chosen by the experimenter and defines the minimum acceptable 
loss of event rate. 




III. THE GRAVITATIONAL WAVEFORM 



The usual restricted PN waveforms are essentially a power series in the velocity v, plus logarithmic terms which 
arise at order and beyond. As we are only considering waveforms up to order in this study we can neglect the 
logarithmic terms and deal with a pure power series. We can create the Fade approximation by simply equating the 
power series on the left hand side to a rational function of the right hand side, i.e. 
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jZa.v'^ - PSi - , (17) 

k=l 

such that n + 1 = N + M + 1. For a more complete discussion of Pade approximation, we refer the reader to 
reference [ p3[ . 

While it is possible to write the Fourier domain waveform in terms of the chirp masses, in order to use Fade 
approximation we need to use a more basic format for the wave. Before progressing, it is necessary to clarify some 
nomenclature. The phase of the waveform is normally constructed using a FN-expansion in the velocity, v. As we 
are dealing with what is in effect a Taylor series, we will call any waveform using FN-approximation a T-approximant 
waveform, while any waveform which uses a Fade approximation, we call a P-approximant waveform. 



A. T-approximant Waveforms. 

Using the Stationary Fhase approximation, we write the Fourier transform of the GW as 

M/) = -4r^/'cxp{z [v'- J]}, 

where the phase, ip, is given by the expression 



^ = 27r/io - (/-o + 2 



dv (vi 
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where w/ = {nmfY^^ is the instanteous velocity at frequency /, viso is the velocity at the last stable circular orbit 
(LSO) and the normalization constant A is given by demanding (^h\lij — 1. In the comparable mass case, the orbital 
energy and GW flux to 0{v^) order are given by [p6 20 



1 - 



9 + 7? 

12 



81 - 5777 + 7]^ 

24 



(20) 



Fiv) 



32772^1 



/1247 35 



44711 
9072 



V 336 

9271 

77 - 

504 ' 



65 
18 



V 



/8191 



V 672 24 



(21) 



The position of the LSO is unknown in the case of two comparable mass objects. However, it has been hypothesized 
that a more accurate description of u/so may be given by [p3| 
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(22) 



which approaches the Schwarzchild test-mass limit in the case of 77 ^ 0. This in turn should provide more accurate 
knowledge of the position of the LSO. From this we can define the frequency at the LSO as fiso = vf^^/irm. 



B. P-approximant Waveforms. 



The groundwork for constructing the Fade approximation of the energy and flux has already been laid down P3[ . 
In this section we will briefly outline the main points. As we are dealing with only sub-diagonal Fade approximants 
we can denote the rational function by P™^^, where as we increase the level of approximation the parameter 
e = or 1. The advantage of the sub-diagonal approximants is purely computational in that we can cast the flux in a 
continued fraction form and we only have one new coefflcient to calculate as we go to each next level of approximation. 
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1. The P-approximant Energy Function 



Starting with Eqn ( |20| ) we can construct the Pade approximant energy function 

1/2 



1, 



where x = and ep^{x) is an energy flux related to the binding energy of the system given by 
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where the coefficients are given by 
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Constructing the Pade approximation gives us the expression 
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Once again, taking the derivative of Eqn (|23| ) with respect to the velocity yields 

vrj dep^{x) 
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where 



dx 



and the Pade coefficients, c^, are given by 
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2. The P-approximant Flux 



As we are only dealing with waveforms up to 2.5-PN order and can neglect all logarithmic terms, we can simplify 
the steps taken for converting the flux into a form compatible with Pade theory. Starting with Eqn (|T|) we can write 
the PN flux as 

5 

FtM = Fn{v)Y.c,v\ (30) 
where F^iv) is the Newtonian flux given by 

F^{v) = ^v'v'°. (31) 
5 

With this representation of the flux, the idea is to create a three step mapping to a Pade approximant flux, Fp{v). 
In order to use Pade theory, our power series needs to contain a linear term in the velocity v. The reason for this is 
that the linear coefficient of the Taylor series appears in the denominator of the Pade coefficients. A Taylor series 
like the PN-flux, which contains a null linear term, results in infinite Pade coefficients. The first step therefore is to 
introduce such a term and create a "factored" flux function by the operation 
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where Vpoie is the velocity at the photon ring, = Cg and = c^— Cf._^/vpoie from fc = 1, ..,n. This choice of using 
Vpoie as a normahzing velocity comes naturally from the fact that we expect a pole at r = 3 M in the flux for particles 
in a Schwarzchild metric. The next step is to create the corresponding P-approximant "factored" flux by 



The final P-approximant flux is then given by 



Lk=0 



(34) 



fpjv). 



(35) 



It has been shown that this representation of the flux has faster convergence properties than the usual PN-flux. 



C. Comparison between T and P-approximant Waveforms. 



The question must be asked as to why we should favor one template model over another. While both templates are 
based on the adiabatic approximation, it has been shown that the P-approximant templates are superior in a number 
of ways |23|j29| ]. The main problem with the T-approximant waveforms is that as we go to higher approximations, 
we do not necessarily obtain a better template. In the test-mass approximation, we know that the v'^ and PN 
approximants diverge in opposite directions from the exact numerical flux. This is reflected when we go to the 
comparable mass case in the values of the overlaps between the 2 and 2.5-PN waveforms and the fiducial "exact" 
waveform. For a 10-10 Mq system, the 2-PN waveform achieves an overlap of 0.894, while the 2.5-PN template has 
an overlap of 0.545. Compare this with the results of the P-approximant waveforms : 0.868 for the 2-PN template 
and 0.979 for the 2.5-PN template, and we can already see the benefit of the P-approximant waveforms. Another 
test can be carried out by looking at the Cauchy convergence between the approximations. If we compare the 2 and 
2.5-PN approximations for both types of waveforms, we get a value of 0.496 between the T-approximant waveforms 
and a value of 0.918 for the P-approximant waveforms. 

The next problem encountered with T-approximant waveforms is that it is possible to achieve overlaps of > 0.9. 
The price that has to be paid, however, is a large bias in the parameter values. For example it has been shown that 
for a [1.4,10] AIq system, it is possible to achieve an overlap of 0.94 with a 2.5-PN template However, this is 
achieved using a template with masses of 0.82 and 20.4 Mq. In comparison, the P-approximant template achieves a 
higher overlap (0.994) with lower biases in the parameter values ([1.37, 10.02] AIq). From this we must conclude that 
the PN waveforms are not faithful templates. 

Taking into account how badly the PN templates perform in the test mass case, we should be worried about their 
performance in the case of two comparable masses. On the other hand, it is due to their faster convergence properties 
and lower biases in parameter estimation, why we believe that the P-approximant templates will give superior results 
in the search for gravitational waves. 



IV. THE NUMERICAL METHOD 



Our primary aim in this study is to calculate the number of templates needed for a GW search using P-approximant 
waveforms. This is closely followed with the aim of developing a numerical code which can be used for any type 
of waveform. As the calculation of the metric tensor depends on the derivatives of the waveforms, it is sometimes 
more convenient to do this numerically. The eliminates the need for calculating the derivatives analytically, which 
for example, in the case of P-approximant waveforms is very messy. In order to calculate the numerator of Eqn (|l5| ) 
there are a number of possible methods. Due to the shapes of the parameter spaces (Fig ^), we chose to use a two-fold 
Monte Carlo method for evaluating the double integral in the same equation. This method is a combination of two 
different Monte Carlo methods called "hit-or-miss" and "crude" methods p^~|3^ . The "hit-or-miss" method is used 
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FIG. 2. The To-r2 and tq-ts parameter spaces for a range of total masses m £ [0.4, 1OO]M0. The range of astrophysically 
interesting sources is contained in both cases within the wedge. 



for calculating the fractional area of the parameter space, 5A. Briefly, this involves enclosing the parameter space 
inside a boundary of known area, and throwing points at random inside the boundary. The fractional area is then 
given by SA ^ A x (n/N), where A is the area of the known boundary, n is the number of points falling inside 
the region of interest and N is the total number of points. If the search manifold was Euclidean we would only 
need to calculate the value ^^g once. For a toy model where all waveforms arc terminated at a constant frequency 
irregardless of mass, then it has been shown that the manifold is approximately Euclidean A suitable co-ordinate 
transformation maps the manifold onto a globally Euclidean space where the metric is constant. However, as soon as 
we include the frequency at the LSO, fiso, in the termination of the waveforms, the search manifold becomes a curved 
Riemannian manifold and a globally Euclidean manifold no longer exist |^0j . This is where the "crude" method comes 
in. We can us it to calculate the mean value of ^Jg on the manifold. By using the chirptimes introduced earlier as 
coordinate basis, we have a choice of two coordinate systems within which to work. These are the tq-T2 and tq-t^ 
systems. What we found, from a numerical point of view, was that the fractional area SA converged quickest in the 
T0-T2 coordinate system. From a computational point of view, this makes it a more desirable set of coordinates as it 
decreases the number of necessary calculations. 

When generating numerical waveforms, the Nyquist theorem states that for a bandwidth limited signal with a critical 
frequency, fNyqi then in reconstructing the time series waveform our sampling rate, fsamp, must be at least twice 
the critical frequency. This is to ensure that any frequency above the critical frequency will not lead to an 'aliasing' 
in the waveform. This is when higher frequency components of the wave occur within the required bandwidth. As 
previously stated, in this investigation it was decided to evolve the waveforms as far as fiso or to some arbitrary cutoff 
frequency. The choice is made as follows. When using the arbitrary cutoff, a frequency of fcut — 2 kHz is chosen. 
This means that all binaries with total masses in the range m G [2.1, 100] Mq are evolved up to and beyond their 
natural cutoff frequencies (by this it is assumed that the natural cutoff frequency for an inspiral waveform is at the 
last stable circular orbit), while only the lowest mass binaries have their waveforms cut off early. It is not expected 
that this will have any great effect on the number of templates as 99% of the waves energy is extracted by the time 
the GW frequency has reached 750 Hz [^. This can easily be seen from looking at the power spectrum of h{f ) which 
is proportional to 1'^ . 

To ensure there is no aliasing in the construction of our waveforms we have chosen the sampling frequency to be 
fsamp =^ 6 kHz if fiso > 1.5 kHz, or, fsamp — ^fiso if fiso < 1.5 kHz. As our waveforms have a maximum bandwidth 
of 2kHz, our sampling frequency is always at least three times fNyq- To cut down on computation time we are 
starting out however with the Fourier transforms rather than the time series waveforms. 

In order to investigate the reliability of the numerical code, we tested it against known results calculated semi- 
analytically. These results were calculated previously for RPNFT waveforms up to 2-PN order using older noise 
curves [ 
minima' 



8|. The basic assumptions of that study were as follows : a rectangular template lattice was used with a 
match of 0.97 and a SNR of 1. It was assumed that the minimal individual mass was 0.2 Mq. Using these 
same assumptions, we calculated the template number for identical waveforms using the numerical method outlined 
above. The results from the numerical code coincided with previous results to within 5%. Thus, we can be confident 
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about our results for template number for both the T and P-approximant waveforms using the new noise curves 



V. RESULTS AND DISCUSSION 



For the initial ground-based intcrfcromcters the most likely sources to be detected are solar-mass BH-BH binary 
systems. This is because these systems have an LSO frequency in or around the most sensitive points of the detector 
noise curves and are therefore the best candidates for detection. For completeness, however, we would also like to 
know how many templates are needed to search for BH-NS and NS-NS binaries as well. With this in mind, we 
calculated the number of templates for two separate cases. In the first case we chose a minimum individual mass of 
1 Mq, while in the second we chose 5 Mq. The maximum total mass in each case is different due to the different lower 
frequency cutoffs of the detectors. The values used are 86 Mq (LIGO, GEO), 100 (VIRGO) and 50 (TAMA). 
For both cases we used the same assumptions. Namely, we choose a minimal match of AiAi = 0.97. This means that 
we have a loss in detection rate loss ~ 10%. All waveforms are terminated at the LSO frequency or 2 kHz, whichever 
is reached first. And finally we choose a rectangular lattice of templates. 

As well as template number, we need to know the computational cost of the storage of these templates. Schutz Q 
has provided a relation between the number of templates and the CPU requirements for a GW search. In brief terms, 
the data streams from the detector are broken into segments of D real numbers. The length of the longest filter is F 
real numbers, such that F >> D. To filter the data through M templates of length F requires N D (16 + 3 log2 F) 
fioating point operations (flops). If the sampling frequency is 2 fjyyq, then the computational power required is 

P~MfNyqm + 6\og.,F) (36) 

flops. If we assume that the length of the longest filter is due to a minimal mass template such that m™*" -I- m™"' = 
m™*", T] = 0.25, then the length of the longest filter, F, is given by 

fNyg [irfion.)-'" (2m™") (37) 

where we once again assume a sampling rate of 2/jvi/g. These equations can be amended accordingly depending on 
the sampling rate used. 

The number of templates needed for a one-step GW search using 2 and 2.5-PN waveforms are presented in Figs § - 
while the computational cost is given in Figs |^ - ^. We present the results for these particular approximations as we 
believe that they currently represent the best templates for a GW search. We can see by comparing between sets of 
results, that in most cases we need more templates to search for inspiralling binaries using P-approximant waveforms 
than we do using T-approximant waveforms. Initially, this may seem counter-intuitive. We would normally assume, 
that if the P-approximant waveforms do indeed bare a closer resemblance to the true waveform of a signal, we would 
need a lower number of templates. We thus need to explain the large increase in the number of templates as we go 
from T to P-approximant waveforms. 

In order to show why we need more templates using P-approximant waveforms, we carried out two experiments. 
In the first experiment we use the fact that the extrinsic components of the metric remain constant over the entire 
manifold. It is only the intrinsic components that change. We therefore held all of the waveform parameters constant 
except for either tq of T2. We chose a number of different points on the search manifold and varied the intrinsic 
parameters by 10~^ seconds. We then calculated the proper distance between both templates. We know from 
Eqn (|l^) that there is a connection between the overlap of two nearby templates and the proper distance between 
them. To first order approximation in AA'^ , the proper distance is given by 



h{\^') /i(A^ -t- AAA')^ (38) 

What we find is that the proper distance between two P-approximant templates is greater than the proper distance 
between two T-approximant templates at constant points on the manifold. This is explained by the fact that not only 
does the proper area of the manifold change as we go to successive approximations, but it also changes as we go from 
one type waveform to another. This can be seen from examining the eigenvalues of the metric tensor as we firstly go 
from one approximation to another, and secondly, by using different types of waveforms. For the second experiment, 
we calculated the span of each template at different points on the manifold. We chose a number of different mass 
ranges on the manifold. In Fig |^ we present the results for the choice of [1, 1] Mq, [1, 10] Mq and [10, 10] Mq while 
changing the Newtonian chirptime tq. The solid curve represents the span of the T-approximant template, while the 
dashed curve represents the span of the P-approximant template. We can see in the case of [1, 1] Mq, there is not 
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much difference in the span of the templates when using either of the approximations. However, as we move to higher 
and higher masses, the P-approximant templates have a smaller span on the search manifold than the T-approximant 
templates. Once again, this can be explained by looking at the eigenvalues of the metric. The orientation and span 
of each template is dependent on the eigenvalues of the locally diagonal metric. We can deduce from this that it is 
in the high mass range of the manifold that is mainly responsible for the change in proper area as we go from T to 
P-approximant waveforms. This fact, combined with the orientation of the template, sufficiently explains the need 
for extra templates when using the P-approximant waveforms. The numerical results are presented in Appendixes B 
and C. 

We should not be discouraged by the vast increase in templates for P-approximant waveforms at certain approxi- 
mations. As we know that the 2.5-PN P-approximant template is vastly superior to its T-approximant counterpart, 
what is most important is the comparison between 2-PN templates. In both cases that we studied, there is only a 
small increase in template number needed as we go from T-approximant to P-approximant waveforms. This increase 
is more than acceptable due to the fact that the P-approximant templates are more robust and reliable in general. 
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FIG. 3. A comparison of the number of templates needed at 2 and 2.5-PN orders for a GW search using T and P-approximant 
waveforms assuming a minimal match of 0.97 and a minimum total mass of 2Mq. 
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FIG. 4. A comparison of the number of templates needed at 2 and 2.5-PN orders for a GW search using T and P-approximant 
waveforms assuming a minimal match of 0.97 and a minimum total mass of 10 Mq. 
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FIG. 5. A comparison of the computational cost at 2 and 2.5-PN orders for a GW search using T and P-approximant 
waveforms assuming a minimal match of 0.97 and a minimum total mass of 2Mq. 
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FIG. 6. A comparison of the computational cost at 2 and 2.5-PN orders for a GW search using T and P-approximant 
waveforms assuming a minimal match of 0.97 and a minimum total mass of 10 Mq. 
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FIG. 7. The span of the T-approximant (sohd) and P-approximant templates (dashed) on the tq axis. From top to bottom, 
the range of masses is [1, 1] Mq, [1, 10] Mq and [10, 10] Mq. The dot-dashed line corresponds to an overlap of 0.97 



VI. CONCLUSION 



In this study we have calculated the number of templates needed for a one-step GW search using T-approximant and 
P-approximant waveforms. We first displayed how the detection of GW can be examined by evoking pre-established 
techniques of differential geometry. By treating the waveforms as vectors in a Hilbert space we can define an N- 
dimensional Riemannian manifold upon which we can define a scalar product, vector norm and a metric tensor. By 
choosing a particular parameterization for the phase of the GW we showed how the search for GW need only take 
place on a 2-D submanifold defined by a set of independent parameters intrinsic to the binary system in question. 
The operation of projecting onto this sub-manifold has a two-fold interpretation. From a geometric point of view it 
is a minimization of distance between two templates, while from a data analysis point of view it is a maximization of 
parameters. We have also outlined the development of an easily modifiable numerical code for calculating template 
number. This code uses the latest analytical models of the noise spectral density of the initial ground-based detectors. 

This code uses a combination of "crude" and "hit-or-miss" Monte Carlo methods to calculate the proper area of 
the search parameter space. We found that there was a particular coordinate system {T0-T2) which seemed to be 
particularly advantageous from a computational point of view. At present the code works for waveforms up to 2.5-PN 
order. We outlined the difference between the T and P-approximant comparable mass waveforms, as well as how to 
construct the P-approximant energy and flux functions. These waveforms were generated to the LSO for waveforms 
with //so < 1.5 kHz. For waveforms with //so > 1.5 kHz we introduced a cutoff frequency of fcut = 2 kHz. This is 
entirely consistent as most of the waves energy has been extracted by ~ IkHz. Although it is BH-BH binaries that 
present the best possibility for detection by the first generation of detectors, we calculated template number for two 
seperate cases assuming minimum individual masses of 1 Mq and 5 Mq 

We found that in general we need more templates when using P-approximant waveforms. In some cases (e.g. 1.5- 
PN order) the number of templates needed increased by a factor of 3. However, the 1.5-PN templates are known to 
be unreliable. The case we were most interested in was the comparison in performance of the 2-PN templates. We 
found that we while we still need a greater number of templates, this increase in small, and when weighed against the 
known robustness of the P-approximant templates, the increase is acceptable. We were able to explain the increase in 
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template number by the fact that the proper area of the parameter space is dependent on the level of approximation 
and type of waveform used. We used this fact to explain how the span of the P-approximant templates decreases 
faster than the corresponding T-approximant templates as we increase the total mass. 

These results are encouraging for the use of P-approximant templates over the usual PN waveforms in the search 
for GW. 
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APPENDIX A: PROPERTIES OF THE METRIC PROJECTION 



We have defined the metric projection as 

7m.=^m.-%^, (Al) 

>0 

What we will now show is that there is a two-fold interpretation of the projection depending on whether one is looking 
at the problem from a geometrical or detection point of view. As one follows on closely from the other, we will first 
examine the meaning of the projection from a geometrical point of view. We firstly define the distance between two 
points (or templates) on a 4 — Z) manifold as 

ds^ = AA'^AA^ (A2) 



= e„„ AA°AA" + 2e„, AA°AA^' + AA^AA^ (A3) 

where we use arbitrary coordinates A^ = {x, y, z, w}. Any function is minimized or maximized if its first derivative is 
zero. Therefore, the proper distance, ds^, is minimized if 

Initially we want to project from an n — D space to an m — _D orthogonal sub-space {m < n). We can project onto the 
orthogonal sub-space as follows. We first choose to minimize the proper distance with respect to the a;-coordinate. 
Therefore, taking the first derivative with respect to AA°, where AA° = jx — xo| 

^ 'ds') = AA" + 2^0, AA^" = 0, (A5) 



9AA0 



The trivial solution is if AA = AA-' = 0. Therefore, Eqn (A5) implies 

^ Lo = Co. = 0. (A6) 
The proper distance between points on the n — 1 sub-manifold is then given by 

= AA^AA^ (A7) 



where jf^i, is the metric tensor on the 3 — D submanifold defined by Eqn ( Al ) and da'^ < ds^ . We can therefore see 
how the projection of the metric can be interpreted as a minimization of distance between points on the manifold. 
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FIG. 8. The minimization of template distance by means of projection onto orthogonal subspaces. 



From a detection point of view we have seen in Section || that the overlap between two templates yields the 
expression 

{u{X0) \u{X^ + AX^)) = 1 - i e^.AA^AA^ (A8) 



The left hand side of Eqn (AS) is maximized if the 2nd term on the right hand side is minimized, i.e. 

^ (e,,.AA^AA'^) = 0. (A9) 



a AAA" 



From Eqn (A4), this process is then equivalent to the minimization of template distance shown above. We can thus 
equate the projection of the metric from one subspace onto another as a maximization over the parameters defining 
each subspace. 



APPENDIX B: TEMPLATE NUMBER FOR A ONE-STEP GW SEARCH. 



Detector 


Ta 


P2 


Ta.s 


P2.5 


LIGO 


9.77 X 10^ 


1.07 X 10* 


5.64 X 10^ 


1.47 X 10" 


GEO 


1.57 X 10" 


1.71 X 10" 


9.13 X 10^ 


2.1 X 10" 


VIRGO 


1.26 X lO'^ 


1.38 X 10^ 


7.67 X 10" 


1.62 X 10^ 


TAMA 


1.34 X 10^ 


2.26 X 10^ 


7.92 X 10^ 


2.27 X 10^ 



TABLE I. The number of templates needed to search for inspiralling binaries using T and P-approximant waveforms at 2 
and 2.5-PN approximations assuming a minimal match of 0.97. The minimum individual mass is 1 Mq and the total mass 
range is m e [2, m"""^] Mq 
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Detector 


Ta 


P2 


^2.5 


P2.5 


LIGO 


51 


75 


26 


92 


GEO 


94 


134 


49 


161 


VIRGO 


808 


1.02 X 10* 


439 


1.15 X 10* 


TAM:V 


1() 


:!l 
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TABLE II. The number of templates needed to search for inspiralling binaries using T and P-approximant waveforms at 
2 and 2.5-PN approximations assuming a minimal match of 0.97. The minimum individual mass is 5 Mq and the total mass 
range is m € [10, m""'"] Mq 

APPENDIX C: CPU REQUIREMENTS FOR A ONE-STEP GW SEARCH. 



Detector 


T2 


P2 


T2.5 


P2.5 


LIGO 


2.12 X 10" 


2.32 X 10"* 


1.22 X 10" 


3.19 X 10" 


GEO 


3.4 X 10^ 


3.71 X lO'' 


1.98 X lO'' 


4.55 X 10® 


VIRGO 


3.08 X 10^° 


3.37 X 10^° 


1.87 X 10^° 


3.96 X 10^° 


TAMA 


2.57 X 10* 


4.33 X 10® 


1.396 X 10* 


4.38 X 10* 



TABLE III. The CPU requirements (in flops) needed to search for inspiralling binaries using T and P-approximant waveforms 
at 2 and 2.5-PN approximations assuming a minimal match of 0.97. The minimum individual mass is 1 Mq and the total mass 
range is m € [2, m"""^] M© 



Detector 


Ta 


P2 


T2.5 


P2.5 


LIGO 


1.31 X 10' 


1.93 X 10' 


6.68 X 10" 


2.36 X 10' 


GEO 


2.42 X lO'' 


3.44 X 10'^ 


1.26 X 10'^ 


4.41 X lO'' 


VIRGO 


2.3 X 10* 


2.91 X 10* 


1.25 X 10* 


3.27 X 10* 


TAMA 


3.71 X 10® 


7.18 X lO'' 


1.62 X 10® 


9.73 X 10® 



TABLE IV. The CPU requirements (in flops) needed to search for inspiralling binaries using T and P-approximant waveforms 
at 2 and 2.5-PN approximations assuming a minimal match of 0.97. The minimum individual mass is 5 Mq and the total mass 
range is m € [10, rrT"''] Mq 
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